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Abstract. Two relativistic X-ray jets have been detected with the Chandra X-ray observatory in 
the black hole X-ray transient XTE Jl 550-564. We report a full analysis of the evolution of the two 
jets with a gamma-ray burst external shock model. A plausible scenario suggests a cavity outside 
the central source and the jets first travelled with constant velocity and then are slowed down by 
the interactions between the jets and the interstellar medium (ISM). The best fitted radius of the 
cavity is ~0.36 pc on the eastern side and ~0.46 pc on the western side, and the densities also show 
asymmetry, of ~0.015 cm -3 on the east to ~0.21 cm -3 on the west. Large scale low density region 
is also found in another microquasar system, H 1743-322. These results are consistent with previous 
suggestions that the environment of microquasars should be rather vacuous, compared to the normal 
Galactic environment. A generic scenario for microquasar jets is proposed, classifying the observed 
jets into three main categories, with different jet morphologies (and sizes) corresponding to different 
scales of vacuous environments surrounding them. 
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INTRODUCTION 

Microquasars are well known miniatures of quasars, with a central black hole (BH), an 
accretion disk and two relativistic jets very similar to those found in the centers of active 
galaxies, only on much smaller scales (Mirabel & Radriguez 1999). 

Since discovered in 1992, radio jets have been observed in several BH binary systems 
and some of them showed apparent superluminal features. In the two well known 
microqusars, GRS 1915+105 (Mirabel & Radriguez 1999) and GRO J1655-40 (Tingay 
et al.1995; Hjellming & Rupen 1995), relativistic jets with actual velocities greater than 
0.9c were observed. In some other systems, small-size "compact jets", e.g. Cyg X-l 
(Stirling et al. 2001), and large scale diffuse emission, e.g. SS433 (Dubner et al. 1998), 
were also detected. 

XTE J 1550-564 was discovered with RXTE in 1998 during its strong X-ray outburst 
on September 7 (Smith 1998). It is believed to be an X-ray binary system at a distance of 
^5.3 kpc, containing a black hole of 10.5±1.0 solar masses and a low mass companion 
star (Orosz et al. 2002). Soon after the discovery of the source, a jet ejection with an 
apparent velocity greater than 2c was reported (Hannikainen et al. 2001). In the period 
between 1998 and 2002, several other outbursts occurred but no similar radio and X-ray 
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FIGURE 1. The smoothed Chandra X-ray 
images of the eight observations of XTE 
J1550-564 and the two jets together. The 
green elliptical regions are source emission 
regions by wavdect. Observation 4 shows 
the good alignment of the two jets and the 
central source. 



flares were detected again in these outbursts (Tomsick et al. 2003). 

With the help of the Chandra satellite, Corbel et al (2002) found two large scale X- 
ray jets lying to the east and the west of the central source, which were also in good 
alignment with the central source. The eastern jet has been detected first in 2000 at a 
projected distance of ~21 arcsec from the central black hole. Two years later, it could 
only be seen marginally in the X-ray image, while a western counterpart became visible 
at ^22 arcsec on the other side. The corresponding radio maps are consistent with the 
X-ray observations (Corbel et al. 2002). 

There are altogether eight 2-dimentional imaging observations of XTE J 1550-564 in 
Chandra archive during June 2000 and October 2003 (henceforth observations 1^8). 
Here we report a full analysis of these X-ray data, together with the kinematic and 
spectral evolution fittings for all these observations. 

OBSERVATIONS OF XTE J1550-564 

The basic information of observations 1~8 is listed in Table 1, including the observation 
ID, date, and the angular separation between the eastern and western jets and the central 
source. The positions are obtained by the Chandra Interactive Analysis of Observations 
(CIAO) routine wavdetect (Freeman et al. 2002). In observations 5 and 6, no X-ray 
source is detected by wavdetect at the position of the eastern jet. However, from the 
smoothed images (Fig.l), a weak source could be recognized in observation 6. We 
thus select the center of the strongest emission region as the position of the jet in that 
observation. We calculate the source centroid for the central source and the X-ray jet 
respectively and for all the five observations, the calculated position changed by less 
than 0.5 7/ . Therefore, an upper limit of 0.5" is set for the error of the jet distance. 

From Table 1 and Fig.l, we could see clearly that an X-ray emission source is detected 
to the east of the central source in the first four observations and another source is 
detected to the west in the last five observations. Calculations also show that these 
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TABLE 1. XTE J155Q-564 Chandra Observations 

Angular Separations (arcsec) Powerlaw Fitting for the western jet 



Num 


ID 


Date 


Eastern Jet 


Western jet 


Photon Index 


Flux (ergs cm 2 s 


_1 ) 


1 


679 


2000 06 09 


21.5±0.5 










2 


1845 


2000 08 21 


22.8±0.5 










3 


1846 


2000 09 11 


23.4±0.5 










4 


3448 


2002 03 11 


28.6±0.5 


22.6±0.5 


1.75±0.11 


(1.9±0.4) x 10" 


-13 


5 


3672 


2002 06 19 




23.2±0.5 


1.71±0.15 


(1.6±0.3) x 10" 


13 


6 


3807 


2002 09 24 


29.2±0.5 


23.4±0.5 


1.94±0.17 


(8.6±1.5) x 10" 


-14 


7 


4368 


2003 01 28 




23.7±0.5 


1.81±0.22 


(5.5 ±1.0) x 10" 


14 


8 


5190 


2003 10 23 




24.5±0.5 


1.97±0.20 


(3.1 ±0.6) x 10" 


-14 



two sources, when presented in a single combined image, are in good alignment with 
the central compact object with an inclination angle of 85.9°±0.3°. By calculating the 
average proper motion, an approximate estimate of deceleration could be seen for both 
jets. 

ENERGY SPECTRUM AND FLUX 

Since the emission from the eastern jet has been studied fully (Corbel et al. 2002; 
Tomsick et al. 2003), we mainly focus our spectral analysis on the western jet. The X- 
ray spectrum in 0.3-8 keV energy band is extracted for each observation of the western 
jet. We use a circular source region with a radius of 4", an annular background region 
with an inner radius of 5" and an outer radius of 15", for each observation. Instrument 
response matrices (rmf) and weighted auxiliary response files (warf) are created using 
CIAO programs mkacisrmf and mkwarf, and then added to the spectra. We re-bin the 
spectra with 10 counts per bin and fit them in Xspec. 

The results of spectra fitting with an absorbed power-law model are also shown in 
table 1 . We use the Cash statistic since it is a better method when counts are low. The 
absorption column density is fixed to the Galactic value in the direction of XTE J 1550- 
564 obtained by the radio observations (Nh = 9 x 10 21 cm -2 ) (Dickey & Lockman 
1990). Our results are quite consistent with previous works by Karret et al.(2003). The 
calculated absorbed energy flux in 0.3-8 keV band is comparable to the value of the 
eastern jet. The observed flux decayed rather quickly, from ^ 1.9 x 10 _13 erg cm -2 s _1 
in March 2002 to only one sixth of this value in October 2003 (see section 4.2). 

JET MODEL 

Kinematic Model 

In the external shock model for afterglows of GRBs, the kinematic and radiation 
evolution could be understood as the interaction between the outburst ejecta and the 
surrounding ISM. Microquasar jet systems are also expected to encounter such interac- 
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tions. In this section, we describe our attempts after Wang et al. (2003) in constructing 
the kinetic and radiation model based on these theories. 

We adopt the model of a collimated conical beam with a half opening angle 0j 
expanding into the ambient medium with the number density n. The initial kinetic 
energy and Lorentz factor of the outflow material are Eq and To , respectively. Shocks 
should arise as the outflow moves on and heat the ISM, and its kinetic energy will turn 
into the internal energy of the medium gradually. Neglect the radiation loss, the energy 
conservation function writes (Huang, Dai, & Lu 1999): 

(r - l)M c 2 + o(Fl - l)m sw c 2 = E (1) 

The first term on the left of the equation represents the kinematic energy of the 
ejecta, where F is the Lorentz factor and Mo is the mass of the original ejecta. The 
second term represent the internal energy of the swept-up ISM, where T sh and m sw 
are the corresponding Lorentz Factor and mass of the shocked ISM respectively, and 
m m = (4/3)nK*m p n(0j/4). 

Coefficient a differs from 6/17 to 0.73 for ultrarelativistic and nonrelativistic jets. 
We adopt the approximation of a ^0.7 after Wang et al.(2003). Equation (1) and the 
relativistic kinematic equations 

dR jB(I> dR jB(I> 

K dt } * l-j8(r)cos0' l ^ ;r l+j8(r)cos0 K) 

can be solved and give the relation between the projected angular separation fi and time 
t. In equation (2), the subscript a and r represent the approaching and receding jets in 
a pair of relativistic jets respectively. R is the distance between the jet and the source, 
which can be transformed into the proper motion separation by \i — R sin 0/5.3 kpc, and 
6 is the jet inclination angle to the line of sight. We can get the \i — t curve numerically 
with the above equations. To be consistent with the work done to the eastern jet, we 
choose the same initial conditions that Tq = 3, = 3.6 x 10 44 erg, and 6j = 1.5. Then 
the parameters needed to be fit are n and G. 

In the case of the eastern jet, the number density of the ISM was assumed as a constant 
in the whole region outside the central source. This assumption does not work well in 
the case of its western counterpart. The western jet decelerated quite fast, requiring a 
local dense environment, but if the ISM is dense everywhere, the jet will be unable to 
travel that far from the central BH. As a result, we consider a model that the ISM density 
varies as the distance changes. For simplicity, we test the ideal case that the jet travelled 
first through a "cavity" with a constant velocity and then through a dense region and was 
decelerated there. A new parameter r, the outer radius of the cavity, is introduced and 
the ISM number density is set to be a constant n outside this region and zero inside. The 
fittings improved a lot but not well constrained because of the limited number of the data 
points. A combination of Lightcurve fitting is required to help the determination. 
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Radiation Model 



In the standard GRB scenario, the afterglow emission is produced by the synchrotron 
radiation or inverse Compton emission of the accelerated electrons in the shock front 
of the jets (Wang et al. 2003 and references there). Wang et al.(2003) found that the 
reverse shock emission, originating from the electrons of the jet when a shock moves 
back through the ejecta, decay rather fast and describe the data of the eastern jet quite 
well. We thus take this model in our work as well. 

Assuming the distribution of the electrons obeys a power-law form, ny e dy e = Ky~ p d% 
for y m < y e < y M , the volume emissivity at frequency v f in the comoving frame is given 
by 

where *i(v,vi,\0 = S^j^F(x)^l 2 dx, with F(x) = xg~K 5/3 (t) and K 5/i (t) is 
the Bessel function. The physical quantities in these equations include q and m e , the 
charge and mass of the electron, B±, the magnetic field strength perpendicular to the 
electron velocity, and v' m and v£, the characteristic frequencies for electrons with y m and 

Assuming the reverse shock heats the ejecta at time to at the radius 7?o, the physi- 
cal quantities in the adiabatically expanding ejecta with radius R will evolve as y m = 

y m (fo)*o/* > 7, = YMRo/R and K = K{t )(R/R )-^ p \B ± = B ± (t )(R/R )- 2 , where 
the initial values of these quantities are free parameters to be fitted in the calculation. 

With these assumptions, we can then calculate the predicted flux evolution of the jets. 
The comoving frequency v' relates to our observer frequency v by v = Dv', where D is 
the Doppler factor and we have D a = 1 /r( 1 — j8 cos 6) and D r = 1 /r( 1 + j8 cos 6) for the 
approaching and receding jets respectively. Considering the geometry of the emission 
region, the observed X-ray flux in 0.3-8 keV band could be estimated by 



(0.3-8 keV) = £\^(*)ARD 3 j v ,]dv, (4) 



where AR is the width of the shock region and is assumed to be AR = 7?/ 10 in the 
calculation. 

To reduce the number of free parameters, we set y m = 100 in our calculation because 
the results are quite insensitive to this value. We choose the time that the reverse shock 
takes place according to our kinematic model in section 3.1. Then we fit the data to find 
out the initial values of K and B±. 

Next step, we combine the kinematic and radiation fitting together. We know that the 
energy and the number density of the gas in the pre-shock and post-shock regions are 
connected by the jump conditions n' — £(T)n and e r — 7] (F)nm p c 2 , where £(r) and 7] (T) 
are coefficients related to the jet velocity. Therefore if we assume the shocked electrons 
and the magnetic field acquire constant fractions (£ e and £ B ) of the total shock energy, we 
have y m = 8 e (/?-2)m p (r- l)/(p- l)m B , K = (p- l)n'yP-\ and£ ± = ^8^7. 

If we further assume that the £ e of the eastern and the western jets is the same, we 
may infer that K oc e f oc n for the two jets. As a result, we search for the combination 
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FIGURE 2. Left: Proper motion fitting with asymmetric cavity. Parameters: 0=68°, r e =14 arcsec, r w =18 
arcsec, n e =0. 00675 cm -3 , ft w =0.21 cm -3 . Right: Reverse shock emission fitting to the X-ray light curve 
of the jets. ^=0.09 cm" 3 , K w =03 cm" 3 , B=0.5 mG, B w =0A mG. 



of parameters that could satisfy the kinematic and radiation fitting, as well as the 
relationship KjK w ~nJn w . 

A set of parameters has finally been found (Please refer to the left panel in Fig. 2). 
The boundary of the cavity lies at r ^14 arcsec to the east and ^18 arcsec to the 
west of the central source. The corresponding number density of the ISM outside this 
boundary is ^0.00675 cm -3 and ^0.21 cm -3 , respectively. These values are both lower 
than the canonical ISM value of ^1 cm -3 , although the value in the western region 
is much higher than in the eastern region. The electron energy fraction relationship is 
satisfied as KjK w ~nJn w ~ 0.03. But the other relation concerning the magnetic field 
strength could not be satisfied simultaneously by these parameters. Although the cavity 
radius and the number density are allowed to vary significantly, the best fitted magnetic 
field strength remains quite stable(~0.4-0.6 mG). One possible interpretation for this 
is that the equipartition parameter varies as the physical conditions of the jet varies; 
an alternative explanation may involve the in situ generation (or amplification) of the 
magnetic field. 

CONCLUSION AND DISCUSSIONS 

External shock model shows that a large scale cavity exists outside XTE J 1550-564. 
This model has also been applied to another X-ray transient H 1743-322. Chandra X-ray 
and ATCA radio observations of this source from 2003 November to 2004 June revealed 
the presence of large-scale (~0.3 pc) jets with velocity v/c ~ 0.8 (Rupen et al. 2004; 
Corbel et al. 2005). Deceleration is also confirmed in this system. The external shock 
model describes the data of this source consistently. A cavity of size ^0.12 pc is likely to 
exist, but not very clear in this case. Even if there is no vacuum cavity, the ISM density is 
found to be extremely low(^ 3 x 10 -4 cm -4 ), compared to the canonical Galactic value. 

These studies led us to the suggestion that in microquasars the interactions between 
the ejecta and the environmental gas play major roles in the jet evolution and the low 
density of the environment is a necessary requirement for the jet to develop to a long 
distance. 
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When putting together all the analyses of microquasar jets, we found that microquasar 
jets can be classified into roughly three groups: small scale moving jets, large scale 
moving jets and large scale jet relics. For the first type, the "small jets", only radio 
emissions are detected. The jets are always relatively close to the central source and 
dissipate very quickly, including GRS 1915+105 (Rodriguez & Mirabel 1999; Miller- 
Jones et al. 2007), GRO J1655-40 (Hjellming & Rupen 1995), and Cyg X-3 (Marti et 
al. 2001). The typical spatial scale is 0^0.05 pc and the time scale is several tenths 
of days. No obvious deceleration is observed before the jets become too faint. For the 
second type, the "large jets", both X-ray and radio detections are obtained, at a place far 
from the central source several years after the outburst. Examples are XTE J 1550-564, 
H1743-322, and GX 339-4 (Gallo et al.2004). The typical jet travelling distance for this 
type is 0.2^0.5 pc from the central engine and deceleration is clearly observed. The last 
type, the "large relics", is a kind of diffuse structures observed in radio, optical and X- 
ray band, often ring or nebula shaped that are not moving at all. In this class, some well 
studied sources, Cygnus X-l (Gallo et al.2005), SS433 (Dubner el al.1998), Circinus 
X-l (Stewart et al. 1993) and GRS1758-258 (Rodriguez et al. 1992) are included. The 
typical scale for this kind is 1^30 pc, an order of magnitude larger than the second type. 
The estimated lifetime often exceeds one million years, indicating that they are related 
to previous outbursts. 

From these properties, it is reasonable to further suggest a consistent picture involving 
all the sources together. We make a conjecture that large scale cavities, exist in all 
microquasar systems. The "small jets" observed right after the ejection are just travelling 
through these cavities. Since there are few or none interactions between the jets and the 
surrounding gas in this region, the jets travel without obvious deceleration. The emission 
mechanism is synchrotron radiation by particles accelerated in the initial outburst. The 
emissions of jets decay very quickly and are not detectable after several tenths of days. 
In some cases (e.g. XTE J1550-564), the cavity has a dense (compared to the cavity) 
boundary at some radius and the interactions between the jets and the boundary gas heat 
the particles again and thus make the jets detectable again. Those are the "large jets". 
The emission mechanism then is synchrotron radiation by the re-heated particles in the 
external shocks. Then, after these interactions, the jets lost most of their kinetic energy 
into the ISM gradually, causing the latter to expand to large scale structures, the "large 
relics", in a comparatively long time (several millions of years). 

The creation of the cavities is not clear at this stage. Possible mechanism may involve 
supernova explosion, companion star winds or disk winds. Since some of the sources 
most likely never had supernovae before and the winds from the companion stars are not 
strong enough, the accretion disk winds may be the most plausible possibility. However, 
these assumptions all require further observations to justify. 

Microquasars are powerful probes of both the central engine and their surrounding 
environment. More studies of the jets behaviors may give us information on the ISM gas 
properties, as well as the ejecta components. It will provide insights of the jet formation 
process and offer another approach into black hole physics and accretion flow dynamics. 
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